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Abstract 

In this work are reviewed several aspects of gravitational lensing produced by astrophysical 
bodies that strongly curve the spacetime in their vicinity. When an object with a photon sphere 
(e.g. a black hole) is interposed between a source and an observer, in addition to primary and 
secondary images due to small deflections of the light rays, there are two infinite sets of images 
called relativistic, produced by light passing close to the photon sphere, which undergoes a large 
deflection. The positions and magnifications of relativistic images can be obtained numerically 
or analytically, in the latter case by an approximate method called strong deflection limit. This 
method is applied to the study of lensing by non-rotating black holes and also by rotating ones. 
For the latter, is also of interest to study the deformation of the shadow cast by the black hole, 
which depends on the rotation parameter. 

PACS number(s): 98.62.Sb, 97.60.Lf, 98.62.Js 
Keywords: gravitational lensing, black hole physics 

1 Introduction 

The presence of a massive body produces the deflection of light passing close to the object, accord- 
ingly to the theory of General Relativity. In 1936, A. Einstein [1] pointed out that the deflection 
of the light of a background star by another star produces two images which can have a great 
magnification in the case of high alignment of the observer, the lens, and the source. He also found 
that the angular separation of the images was too small to be resolved by the optical telescopes 
available at that time. The discover of quasars in 1963 opened up the possibility of really observing 
gravitational lensing effects. These objects, situated at cosmological distances, are very bright and 
have a compact optical emitting region. The magnification can be large and the images are well 
separated in some particular cases. The first gravitational lens discovered, by Walsh et al. in 1979, 
was the quasar QSO 0957+561 A,B. The weak deflection theory of gravitational lensing, devel- 
oped, among others, by Y. G. Klimov, S. Liebes, S. Refsdal, R. R. Bourossa, and R. Kantowski, 
has been successful in explaining the astronomical observations up to now. This theory is based on 
a first order expansion of the small deflection angle. For a detailed treatment see [2], and references 
therein. 

The theoretical research of gravitational lensing by compact objects with a photon sphere, 
such as black holes and naked singularities, has received great attention in the last decade, mainly 
because of the strong evidence about the presence of supermassive black holes at the center of 
galaxies, including the Milky Way [3]. For these lenses, large deflection angles are possible for 
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photons passing close to the photon sphere. These photons could even make one or more complete 
turns, in both directions of rotation, around the black hole before eventually reaching an observer. 
As a consequence, two infinite sequences of images, called relativistic, are formed at each side of 
the black hole. An useful analytical method for obtaining the positions, magnifications, and time 
delays of the relativistic images corresponding to black holes as gravitational lenses, is the strong 
deflection limit, a logarithmic approximation of the deflection angle for light rays passing close to 
the photon sphere. Introduced by Darwin [1] for the Schwarzschild geometry, it was rediscovered 
several times [SHE]; extended to the Reissner-Nordstrom spacetime [9], and to any spherically 
symmetric black holes [10]. Numerical studies of black hole lenses were done too |11H13| . Black 
holes with spherical symmetry coming from alternative theories |14pi5| , string theory or braneworld 
cosmologies |16H24| and even naked singularities [25,26] were considered as lenses. Retrolenses were 
also studied |27H31j . Particular interest have received the supermassive Galactic black hole as a 
possible lens, with the stars in its neighborhood as light sources |31H34j . For other related works 
see Refs. [35rl40| . Kerr black hole lenses were analyzed by several authors [41H47j . Rotating black 
holes present apparent shapes or shadows with an optical deformation due to the spin [4"8rf50] . 
instead of being circles as in the case of non-rotating ones. This subject have been re-examined in 
the last few years [4l ^l5l"rf59] . with the expectation that the direct observation of black holes will be 
possible in the near future [54 ^ 155 1 160] . so the study of the shadows will be useful for measuring the 
properties of astrophysical black holes. Optical properties of rotating braneworld black holes were 
studied in Refs. [61,62j. Other interesting topics are discussed in the recent review articles |63p64j . 

In this paper, several aspects of strong deflection gravitational lensing due to black holes or 
naked singularities are reviewed. The first part of this work corresponds to lenses with spherical 
symmetry and the second part to rotating lenses. In the latter case, particular attention is given to 
the deformation of the shadow due to the spin. Finally, a brief discussion about the observational 
prospects is given and a summary is made. Units such as G = c = 1 are adopted. 

2 Spherically symmetric lenses 

The lensing scenario to be analyzed here consists of a non-rotating object with a photon sphere, 
a point source of light and an observer. The geometry is supposed asymptotically flat, with the 
source and the observer both situated far from the object, in the flat regiorQ. 

2.1 Deflection angle 

For a general spherically symmetric and asymptotically flat metric, having in Schwarzchild coordi- 
nates the form 

ds 2 = -f(r)dt 2 + g{r)dr 2 + h{r){d6 2 + sin 2 9dip 2 ), (1) 
the deflection angle a is given by [25 , 65J 
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where ro is the closest approach distance. The photon sphere radius r ps , corresponding to the 
unstable circular photon orbits, is the largest positive solution 35, 3(> of the equation: 
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Figure 1: Schematic lensing diagrams, where the point light source (S), the deflector (L), the 
observer (O), and the image (I) are shown. The angular diameter distances are <ios (observer- 
source) , cf ol (observer- lens) , and c?ls (lens-source) ; f3 and 9 are respectively the angular source and 
image positions. From left to right, the first two diagrams (a) and (b), correspond to standard 
lensing, with the lens between the observer and the source, and the other two (c) and (d), to 
retrolensing. 

for Schwarzschild geometry r ps = 3M. 

There are two cases where the deflection angle a can be approximated by simple expressions: 

• Weak deflection limit: if vq 3> r ps , the deflection angle can be approximated by the first non 
null term of the Taylor expansion in 1/ro (e.g. Ref. [2j). This limit is suitable for the analytic 
calculations corresponding to the primary and secondary images, formed by photons passing 
at large distances from the deflector. 

• Strong deflection limit: the deflection angle a(ro) diverges when r$ — > r ps . For < ro/r ps — 
1 1, it can be approximated by a logarithmic function [10J: 



where a\ and a>2 are constants that depend on the specific form of the metric functions (see 
Appendix); for the Schwarzschild black hole the values can be calculated exactly [8]: a\ = 2 



to the two sets of strong deflection (or relativistic) images, formed by photons passing close 
to the photon sphere. The strong field limit provides an excellent approximation for the 
deflection angle (see discussion in Refs. [H1I37]). 

In Fig. [H different possible lensing situations are shown. The optical axis is defined as the line 
joining the observer and the lens. The angular positions of the source and the images, seen from 
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the observer, are (3 (taken positive) and 9, respectively. The observer-source (dos), observer-lens 
(doh) and the lens-source (cZls) distances are measured along the optical axis and they are taken 
much larger than the horizon radius. From left to right, the first two plots (a) and (b), correspond 
standard lensing, i.e. the deflector between the source and the observer, while the other two (c) 
and (d), represent retrolensing scenarios |27,29j. The first diagram (a) corresponds to photons 
with ro 3> r ps , which form the primary and secondary images (only one of them is shown). The 
second plot (b) shows a light ray which makes a complete turn around the lens before reaching 
the observer; in this way the first relativistic image is formed ( there is another one due to light 
rays passing by the other side of the lens, not shown). Other relativistic images are obtained with 
photons making more turns around the deflector, not shown in the diagram. The third (c) and 
four (d) plots show one retrolensing image produced by photons with a deflection angle close to it 
(light rays passing by the other side of the lens form another image, not shown). More retrolensing 
images are formed by photons that make complete turns around the lens, with a deflection angle 
close to an odd multiple of tt, not shown in the diagrams. 



2.2 Primary and secondary images 

In the case of standard lensing with a very small source angle /?, the primary and secondary images 
are present, for which the deflection angle a is very small. Then, the angular positions of the images 
are also small. 



Image positions 

It is not difficult to see from the lens diagram, that the lens equation for small lensing angles takes 
the form [2]: 

p = o-!h* a . (5) 

"OS 

The first step is to calculate the deflection angle a in the weak deflection limit, by Taylor expanding 
the integrand of Eq. ([2]) and performing the integral. To first order in 1/ro, one finds that a ~ x/ r 0; 
where x ls a constant which depends on the metric considered; for Schwarzschild geometry x = 4M. 
Using that ro = doi^O and inverting the lens equation, the angular position of the images are: 



\(p±y/p + 4tf%), (6) 



with 



the Einstein radius. For perfect alignment, instead of two point images, an Einstein ring with 
angular radius Oe is obtained. 



Magnifications 

Gravitational lensing conserves surface brightness (e.g. Ref. [2j), then the magnification is the 
quotient of the solid angles subtended by the image and the source 
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which, using Eq. ©, gives for each image: 
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and the total magnification is 



2(3^(3^ + 49% 
(3 2 + 29 2 E 



2 2 



(9) 



(10) 



The magnifications are proportional to 1/(3, so the magnifications grow as the objects are more 
aligned. 

Time delay 

Photons that form the primary and secondary images follow different paths and the time delay 
between them is [21 
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where Zd is the redshift of the deflector. When (3 = an Einstein ring is obtained, so there is no time 
delay for perfect alignment. Large time delays can be obtained if /3/#e 3> L but the magnification 
of the primary image is close to one and the secondary image is very faint in this case. The time 
delay between the images can be a useful observable quantity in the case of transient sources. The 
optimal situation for a variable source is when (3/6~e is small enough to have large magnifications 
of both images, but not too close to zero, so the time delay can be longer than the typical time 
scale of the source. 

2.3 Relativistic images 

Light rays passing close to the photon sphere will suffer a large deflection angle. As a consequence, 
the relativistic images appear, produced by photons making one or more turns around the deflector. 
The lensing effects are more important when the objects are highly aligned, so here is studied in 
detail only this configuration, which is the simplest and more interesting one^|. In this case, the 
angles (3 and 9 are small, and a is close to a multiple of ir. The strong deflection limit is adopted 
for the (approximate) analytic calculations, in which the deflection angle is given by Eq. @, with 
the coefficients a\ and 02 depending on the specific form of the metric. 



Image positions 

The impact parameter b is more easily connected to the lensing parameters than the closest approach 
distance tq. They are related by the equation [251 165] 
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By Taylor expanding Eq. ([12]) and using Eq. © to obtain that b'(r ps ) = 0, it is straightforward to 
see that 

b - b ps k(r - r ps ) 2 , (13) 



2 For a more general treatment, without the small angle approximation, see Ref. [31| . 
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with bp S = [h{r ps )/ f(r ps )\ l l 2 and k = b"(r ps )/2 (for Schwarzschild lens, b ps = 3\/3M, k = 
\/3/(2M)). Then, by replacing Eq. <H3J> in Eq. @J, one obtains 

a{9) = -a In ( A - + c 2 + O f A - lV (14) 



where 

ci = | (15) 



and 



For Schwarzschild metric, the values are [TO] : a = 1 and c 2 = ln[216(7 - 4\/3)] - vr. Photons 
having an impact parameter b slightly larger than the critical value b ps will spiral out, eventually 
reaching an observer after one or more turns around the black hole, and those with b smaller than 
b ps will spiral into the black hole, not reaching any observer outside the photon sphere. From the 
lens geometry, it can be seen that 

b = d L sin 9 as d OL 9; (17) 
so, using Eqs. (|17p and (|18p . the strong deflection limit takes the form 



a(9) = - Cl In — - 1 + c 2 + O — - 1 , (18) 

\tlps / \tlps J 

with 8 P s = bp S /doL. The lens equation that relates the lensing angles with the distances can be 
written in the form [llffi 

tan f3 = tan 6 — C3 [tan(a — 9) + tan 9] , (19) 

where C3 = dLs/^OS [H] for standard lensing (Fig. Q](b)) and C3 = dos/doh or C3 = dos/^LS fo r 
retrolensing (Fig. UJ (c) and (d), respectively) [22] • A non- negative angle j3 can be taken without 
losing generality. When f3 7^ two infinite sets of point relativistic images are obtained. [4"]ll lj. For 
high alignment, the deflection angle can be approximated by 

a = ±mir ± Aa m , (20) 

with m£fJ and < Aa m <C 1; m is called the winding number, it is even for standard lensing and 
odd for retrolensing. The plus and minus signs correspond to light rays passing by different sides 
of the black hole. Replacing a in Eq. (|19p . and using that /3 and 9 are small for high alignment, 
the lens equation adopts the form: 

P = m T c 3 Aa m . (21) 



Inverting Eq. (fTHj) to obtain 9(a) 
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and making a first order Taylor expansion around a = rmr, the angular position of the m-th image 
is given by 

9 m = 9° m T( m &a m , (23) 



3 For a throughout discussion about possible lens equations see Ref. [38] 
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Prom Eq. (f2T|) . one has that Aa m = ±(9 m — 0)/c$, and replacing it in Eq. ([23]) leads to 
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which can be written in the form 
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Finally, using that < ( m /c3 *C 1 and keeping only the first order term in Cm/c3, the angular 
positions of the images are given by 



dm — ±# m + (/3 =F ^m) • 

For perfect alignment an infinite sequence of Einstein rings with angular radius 
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is obtained instead of the two sequences of point images. 
Magnifications 

As gravitational lensing conserves surface brightness (e.g. Ref. [2]), the magnification is the quotient 
of the solid angles subtended by the m-th image and the source 
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and using Eq. ([28]) it gives 
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so, to first order in Cn/c3, one has 
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Mm = o^- : -> (32) 

for both sets of images. The first relativistic image is the brightest one, and the magnifications 
decrease exponentially with m. The total magnification for a point source is obtained by adding 
the magnifications of both sets of relativistic images 
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with 5 = for standard lensing and 5 = it for retrolensing. The magnifications of the strong 
deflection images are larger for retrolensing. In both cases, the magnifications are proportional to 
l/c£, which is a very small factor. Then, the relativistic images are very faint, unless P has values 
close to zero, i.e. nearly perfect alignment. For p = 0, the amplification becomes infinite, and the 
point source approximation breaks down, so an extended source analysis is necessary. 
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Extended source 



When the source is extended, the magnification of the images can be obtained by integration over 
its luminosity profile: 

where X is the surface intensity distribution of the source and jl is the magnification corresponding 
to each point of the source. If the source is an uniform disc D(f3 c ,/3 S ), with angular radius j3 s and 
centered in j3 c (taken positive), Eq. (|34p can be written in the form 
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Then, using Eq. (|31|) . the magnification of the relativistic m-th image (with m even for standard 
lensing and odd for retrolensing) when the source is an uniform disc is given by [8] 
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with K(k) = Jq /2 (1 - k 2 sin 2 ^^dcp and E(k) = f* /2 (l - k 2 sin 2 cp) 1/2 dcp the complete elliptic 
integrals of the the first and the second kind, respectively. The total magnification for an uniform 
disc is 

J g e (c 2 +<5)/ci |"j _|_ e (c 2 +<5)/ci _|_ g2vr/cij 

^^Id 2 ^ c^M - 1 ' (38) 

with 5 = for standard lensing and 5 = n for retrolensing. These expressions always give finite 
magnifications, even in the case of complete alignment. 



Time delays 

As pointed out above, the time delay between images is of particular interest in the case of transient 
sources. For the standard lensing configuration, the time delay between the relativistic images [37] . 
resulting from the different paths followed by the photons that form them, at the same side of the 
lens with winding numbers m and n, is given b}Q 
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where k = I^J g(r ps ) / \/ kb ps f(r ps ) , and the plus (minus) sign if both images are on the same 
(opposite) side of the source . For the images at the opposite side of the lens: 
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4 The expressions from Ref. [37] have been rewritten here, by adding the cosmological factor 1 + Zd and expanding 
them to first order in the source position angle, which is measured from the observer instead of from the source, as 



done in Ref. |40] for Schwarzschild geometry. 
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where the image with winding number n is on the same side of the source and the one with m on 
the opposite side. The first term in Eqs. (|39p and (|40p is by large the most important one. The 
time delays between the relativistic images are longer than the time delay between the primary 
and the secondary images. 

3 Rotating lenses 

In the case of rotating black holes, photons follow complicated paths in the vicinity of these compact 
objects. They do not move in a single plane, so the problem is more complex than in the non- 
rotating case. The Hamilton-Jacobi equation determines the null geodesies for a given geometry. 
In the case of Kerr metric the problem can be separated by using Carter constants (e.g. Ref. [49J). 
In other spacetimes this separation can be also done. The caustic structure is also more complex 
than in the non-rotating case [H] . The analytical treatment of rotating black hole lenses is much 
difficult than non-rotating ones. The strong deflection limit was extended to Kerr black hole lenses 
to obtain the positions and magnifications of the relativistic images, to second order in the rotation 
parameter a. For a detailed treatment, see [32] and references therein. 



When a black hole is in front of a luminous background, the light reaches the observer after being 
deflected by the black hole gravitational field; but part of the photons emitted by the source, those 
with small impact parameters, end up falling into the black hole, not reaching the observer, and 
producing a completely dark zone called the shadow. The apparent shape of a black hole is thus 
defined by the boundary of the shadow. For rotating black holes, prograde photons can get closer 
to the black hole than retrograde ones. The shadows present a deformation that grows with the 
rotation parameter, instead of being circles as in the non-rotating case. A schematic example is 
shown in Fig. [2] (left). The deformation also grows with the inclination angle Oi of the observer 
with respect to the rotation axis of the black hole. For a polar observer (9i = or 8i = ir) there 
is no deformation, while for one in the equatorial plane (6i = 7r/2) the gravitational effects on the 
shadow are larger. If there are other parameters associated to the black hole, e.g. the charge, it 
will depend on them too, in a complex way. The contour of the shadow for a given geometry can 
be obtained by following the procedure shown in Ref. [19] for Kerr metric. The first step is to 
separate the Hamilton-Jacobi action for the null geodesies into a radial and an angular part. The 
null geodesies are then parametrized in terms of the conserved quantities £ = L z /E and 7] = Q/E 2 
(with E the energy, L z the axial component of the angular momentum, and Q the Carter constant). 
The quantities £ and r] can be related with the celestial coordinates a and /3, which are the apparent 
perpendicular distances of the image, as seen from the axis of symmetry and from its projection 
on the equatorial plane, respectively (see e.g. @J]). The boundary of the shadow is determined by 
the geodesies of photons with parameters ^ sp h and eta sp h, corresponding to the unstable spherical 
orbits around the deflector. This method was applied to Kerr-Newman |52j and to from alternative 
theories black holes [5T? j loU l lfi2"] . 

For the characterization of the form of the shadow, are useful the observables defined in Ref. |58J: 
the radius R s and the distortion parameter 5 S . The quantity R s is the radius of a reference circle 
passing by three points: the top position of the shadow, the bottom position of the shadow, and 
the point corresponding to the unstable retrograde circular orbit when seen from an observer on 
the equatorial plane. The distortion parameter is defined by D/R s , where D is the difference 
between the endpoints of the circle and of the shadow, both of them at the side corresponding 



Shadows 
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Figure 2: Shadows of a rotating black hole (left) and a rotating naked singularity (right). 



to the prograde circular orbit. The radius R s basically gives the approximate size of the shadow, 
while S s measures its deformation with respect to the reference circle (see [58J for more details). 

Another related topic is the apparent shape of a rotating naked singularity (if these objects 
exist in nature); in this case an open dark arc and a deformed dark disc appear (see e.g. [52|.58j). 
as shown in Fig. (right). The unstable spherical photon orbits with a positive radius result 
in an open dark arc instead of a closed curve; the photons near both sides of the arc can reach 
the observer due to the nonexistence of the horizon. The unstable spherical photon orbits with a 
negative radius constructs the small dark disc: the observer will never see the light rays from such 
directions because they escape into the other infinity by passing through the inside of a singular 
ring. When the observer is on the equatorial plane (9i = vr/2), the same arc exists but the dark 
disc disappears, because the light rays in the direction of negative radius will always hit on the ring 
singularity. Those null geodesies result a straight line. From the arc, two observables can be defined 
citehioki09: the radius R a of a circumference that approximates the dark arc, and the central angle 
(f a , determined by the angle subtended by the arc, seen from the center of the circumference used 
to define R a . 



4 Observational prospects 

Forthcoming instruments are expected to look into the region close to Sgr A*, the supermassive 
black hole at the Galactic center, with mass 4.3 x 1O 6 M and at a distance of 8.3 kpc [3]. The 
apparent angular radius of this supermassive black hole is about 27/ias, so resolutions of the order 
of //as are required for the observation of the lensing effects considered here. More massive black 
holes can be found in other galaxies, but the distances are larger too. In the near future some 
observational facilities, most of them space-based, will be fully operational, and will be able to 
measure in the radio and X bands. Two of them are worth mentioning: 

• RADIOASTRON: space-based radio telescope, launched in July 2011. It will be capable of 
carrying out measurements with high angular resolution, about 1 — 10/ias [54, 55, 60, 66J. 

• MAXIM project: is a space-based X-ray interferometer with an expected angular resolution 
of about 0.1 /xas (see [67] for further details). 

The observation of black hole apparent shapes is a major goal in observational astrophysics, since 
those shadows correspond to a full description of the near horizon region, without any theoretical 
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assumption concerning the underlying theory or astrophysical processes in the black hole surround- 
ings. These instruments will be able to resolve the shadow of the supermassive Galactic black hole, 
which together with other observations, would serve to obtain its parameters in the near future. 
Gravitational lensing would be then an useful tool for the study of black holes. More subtle effects, 
like the comparison with observations of different models of black holes, corresponding to alterna- 
tive theories of gravity, would require a future second generation of instruments. Other interesting 
observational aspects are discussed in Refs. [601153]. 



5 Summary 

The observation of the strong deflection effects will be a very difficult but not impossible task not so 
far in the future. Different topics of strong deflection gravitational lensing produced by a black hole 
or a naked singularity were reviewed in this paper. The weak and the strong deflection limits were 
used for the analytical computation of the positions, magnifications and time delays of the primary, 
secondary and relativistic images, for high alignment, which is the simplest and more interesting 
case. The relativistic images have an angular position about the size of the angle subtended by 
the event horizon and they are strongly demagnified, making their detection a challenge for future 
instruments. The time delay between images for transient sources was also discussed in this work. 
Other relevant aspect is the deformation of the shadow for rotating lenses, which grows with the 
spin, and it would give useful information in the near future, when direct imaging of the Galactic 
supermassive black hole will be possible. 
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Appendix 

In this Appendix, the strong deflection limit coefficients are calculated. The deflection angle can 
be rewritten in the form a(ro) = I(fo) — tt, where I(tq) is the integral in Eq. ([2]). The integral 
I(ro) grows as ro approaches to r ps , where it diverges. Following Ref. [TO], it is useful to define a 
new variable z = [f(r) — /(^o)] / [1 — /( r o)]j an d the functions 



R{Z > ro) = f'(r)h(r) [1 " /(r ° )] (41) 
w(z,r )= - 1 (42) 

v7M - [(i - 7M) z + /(r )] WmrJF 1 

where r = f~ 1 [(l — f{ro))z + f(ro)]. By performing a Taylor expansion of the function inside the 
square root in Eq. (|4"2j) one obtains 

w (z,r ) = 1 ===, (43) 

V V(ro)z + j{r )z 2 

where 

tiro) = ^"(^ [f(ro)h'(r ) - f(r )h(r )] , (44) 
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and 



7(ro) = 2 [fVo)]Wo)] 2 W^ h ^ h '^ " f(ro)f"(r )h(r )h'(r ) 

+/(ro)/'(r ) [^(ro)fc"(r ) - 2[/ l '(r )] 2 ] } . (45) 
With these definitions, the integral /(ro) can be separated into two parts 

I(r ) = I D (r ) + I R (r ), (46) 

where ^ 

Io(ro) = / R(0,r ps )w (z,r )dz, (47) 



o 



and 

»1 



lR(ro) = [R(z,r )w(z,r ) - R(0,r ps )w (z,r )]dz. (48) 
J o 

If ro 7^ r ps we have that 99 7^ and wo ~ l/v^> so the integral Id(j"q) converges. Instead, when 
?~o = r ps) by using Eq. ([3]) one has that ip = and wo ~ l/^i so /d(?"o) nas a logarithmic divergence. 
Therefore, Id(to) is the term containing the divergence at ro = r ps and In(ro) is regular because 
it has the divergence subtracted. By performing the integral in Eq. (|47p and making a Taylor 
expansion, it can be seen that |10| 



where 



I D (r ) = -a 1 ]n( ^-1] + a D + O ( ^ - 1 ) , (49) 



«1 = 7=7 r ' ( 50 ) 

and 

gfeg) , 2[l-/(r pB )] 

«D = ^==ln — r-. 51 

Finally, by Taylor expanding the regular integral Ir{tq) = Ir{t V s) + O(ro/r ps — 1), and defining 
°R = lR,( r ps), the strong deflection limit takes the form shown in Eq. with a\ given by Eq. (|50p 
and 

02 = or + od — 7T. (52) 

In some cases, e.g. in Reissner-Nordstrom spacetime, the integral Ir{t ps ) cannot be calculated 
exactly and it is necessary a numerical integration or a Taylor expansion in terms of a relevant 
parameter to obtain ajj. For an alternative numerical computation of the strong deflection limit 
coefficients see Ref. [9]. 
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